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Abstract 

The effects of line-tying on resistive tearing instability in slab geometry is studied within the 
framework of reduced magnetohydrodynamics (RMHD). It is found that line-tying has a 

stabilizing effect. The tearing mode is stabilized when the system length L is shorter than a 
critical length Lc, which is independent of the resistivity r]. When L is not too much longer than 
Lc, the growthrate 7 is proportional to rj . When L is sufficiently long, the tearing mode scaling 
^ 7^3/5 jg recovered. The transition from 7 ~ to 7 ~ r]^/^ occurs at a transition length 
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I. INTRODUCTION 



The study of the effects of hne-tying on magnetohydro dynamic (MHD) instabihties has 
a long history. In the solar corona, the magnetic field lines of solar coronal loops are deeply 
anchored (line-tied) to the much denser photosphere. It has long been suggested that because 
of the stabilizing effect of line-tying, the loops can remain stable for a much longer period 
of time than expected from MHD instability growth times without line-tying. Raadu was 
the first to analyze the line-tying stabilization of ideal internal kink mode using the energy 



principle. [3| Later, more detailed analyses have been done for various resistive 



and ideal MHD instabi 
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^In recent years, there is a reviving theoretical interest in 
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17| partly due to progress from plasma laboratory 



In this paper, we study the effects of line-tying on resistive tearing instability in one of the 
simplest settings, the slab geometry. The tearing instability in slab geometry is discussed in 
many textbooks. 20, l21[ The equilibrium is a simple current layer where the current is along 
the z direction, and the current density only depends on x. Associated with the current 
layer is a magnetic field By which changes direction across the current layer. This configura- 
tion is ideally stable. With resistivity, the magnetic field can reconnect through the resistive 
tearing instability. To simplify the analysis, usually a strong guide field perpendicular to the 
reconnecting field is assumed such that the incompressible approximation may be justified. 
The configuration we consider here is to place two perfectly conducting end plates perpen- 
dicular to the guide field. Therefore, it is the guide field that is line-tied to the end plates. 
Although this configuration is rather simple, to our knowledge, it has not been studied pre- 
viously. Velli and Hood have studied the line-tied tearing mode in slab geometry; but in 
their case, the conducting end plates are perpendicular to the reconnecting field. The model 
we study here may be regarded as a prototype of more complicated resistive instabilities 
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that arise in coronal heating models based on Parker's scenario. [22, 

The rest of the paper is organized as follows. Sec. [Ill sets up the physical model and 
the governing equations of the system. Sec. IIIII discusses some general properties that can 
be deduced from the equations without actually solving them. Sec. [IVl gives a concise 
description of the numerical methods. The numerical results are presented in Sec. [Vl We 
summarize and conclude in Sec. 
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II. PHYSICAL MODEL AND GOVERNING EQUATIONS 



To simplify the problem, we make the following assumptions: (1) the plasma density is 
a constant; (2) a strong, constant guide field along the z axis is present; and (3) the length 
scales along the guide field are much longer than the length scales in the perpendicular 
directions. Under these assumptions, the system is may be described by the well known 
reduced MHD (RMHD) equations, [l, 2] In this framework, after proper normalization, the 
magnetic field and the velocity can be expressed in terms of the flux function and the 
stream function as B = z + V^?/' x z and u = V_l0 x z. The governing equations are 



dSl + [0, ^\ = d,J + [tp, J] + i^vln, 



(1) 

(2) 



Where v and rj are the viscosity and the resistivity; = + ydy is the perpendicular 
gradient; = — V5_0 , J = —V^tp are the vorticity and the electric current; and (p] = 
dyipdx4> — dxipdycj) is the Poisson bracket. This set of RMHD equations is widely used in 



studying t 
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le coronal heating and current sheet formation problems. 
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34l | Two perfectly conducting end plates are placed aX z = ±L/2 to provide 



the line-tied boundary condition. Although the line-tied condition is an idealization in the 
sense that the solar photosphere is neither completely rigid nor perfectly conducting, it is a 
common approximation in solar physics. Strictly speaking, the RMHD ordering ^ V_l is 
incompatible with the line-tied boundary condition, due to the presence of a boundary layer 
near the end plate where ~ V_l. However, Scheper and Hassam have shown through a 
boundary layer analysis that RMHD can still be used in the coronal volume, as the boundary 
layer only appears when the higher order correction to is calculated. The higher order 
correction to -B^, however, is completely decoupled from ip and 

We assume that the equilibrium only depends on x; i.e., B 
Consider perturbations of the form (p = 
Linearizing the RMHD equations gives 



35| 

z - dx^py = z + By-y. 
X, z) exp{ikyy+'yt), and pj = p){x, z) exp{ikyy+'yt). 



-fV(p = (d, + ikyBy)Vp) - ikyB'yp) + z/P^ 
'^p) = (9^ + ikyBy)(p + if]'Dp). 



(3) 
(4) 



3 



where V = — ky and the prime denotes dx- The hne-tied boundary condition in the 
RMHD approximation is simply (p = a.t z = ±L/2. In this study we will focus on the 
inertial tearing mode; hereafter the viscosity u = is assumed throughout the rest of the 
paper. 

III. ANALYTICAL CONSIDERATIONS 

Although the linearized equations are difficult to solve analytically, some general proper- 
ties can be readily deduced from analytical considerations. First of all, let us assume that 
there exists a critical length Lc such that 7 = 0. By setting 7 = and rescaling z with 
respect to Lc (i.e., z — >• Lcz), Eqs. ([3]) and (jlj) become 

- -^dz(f) = ikyBycj) + ?7Pi/', (6) 

with the boundary condition 0|z=±i/2 = 0. Eqs. ([5]) and ([6]) may be regarded as an eigen- 
value problem with l/L^ as the eigenvalue. We further notice that by letting rjcf), the 
parameter 77 completely drops out of the equations. Therefore, we conclude that the critical 
length Lc is independent of rj. Some caveats should be mentioned: (1) the existence of is 
not guaranteed; (2) the marginal stability may not occur at 7 = 0, as pure oscillatory modes 
are not precluded by the equations. Bearing these caveats in mind, we found numerically 
that Lc does exist, and indeed it is independent of rj. 

Now let us consider an unstable system not too far from marginality, i.e., L > Lc- If the 
unstable mode is so slowly growing, we may neglect the inertia term (the left hand side) in 
Eq. ([3]). Physically, that means that the system remains force- free throughout the whole 
domain. Neglecting the inertia term in Eq. ([3]), and dividing Eq. (jlj) by 77, we come to the 
following equations for the force-free approximation: 

= {d,+ lkyBy)Vi, - ikyBlij, (7) 

(7/r/)^A = (a, + tkyBy) (0/r/) + V^. (8) 

It immediately follows that in the force-free approximation, both the growthrate 7 and the 
stream function scale as rj. Since the inertia term ~ 70 ~ 77^, it is expected that the 
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force-free approximation works better for smaller 77, when everything else remains the same. 
Recall that in a periodic system, the magnetic field is approximately force-free only in the 
"outer" region, i.e., far from the resonant surface where k ■ B = 0. In the inner layer, the 
resistivity and the inertia are both taken into account. The inner layer thickness 5 ~ r/^/^, 
and the growthrate 7 ~ rj^^^ as a result. {sg], [s^I Naturally, we expect the 7 ~ 77'^/^ scaling 
to be recovered if the system is sufficiently long. The operative question is, at what length 
does this transition occur? We will address this question later. 

Finally, we would like to say a few words about the boundary condition. Both Eq. ([3]) 
and Eq. (jlj) are first order differential equations in z, for ip and respectively. Therefore 
two boundary conditions are needed. We have two boundary conditions (f)\z=±L/2 = for 0, 
and none for ip. If -0 is given, then one may integrate Eq. (jl]) along z to solve (p. In general, 
the (j) so obtained can not satisfy both boundary conditions, unless some solubility condition 
for ijj is satisfied. Acting on (j4]) by J^^^^dz e^'''"^"^ , after integrating by parts and using the 
boundary condition on 0, we have the solubility condition 

pL/2 ^ pL/2 

7/ dze'^^^«'^ = 7] dze'^y^y'V^j. (9) 

J~L/2 J-L/2 

This solubility condition is not used in the rest of the paper, but we have used it to check 
the accuracy of the numerical solution. 



IV. NUMERICAL METHOD 



The numerical method we employ may be regarded as a variation of the method proposed 
by Evstatiev et al. in Ref. l^]- In essence, this method is a shooting method in two 
dimensions. We may rewrite Eqs. ([3]) (HI) in the following form as an evolution equation 
along z: 

d_ 

dz 

where the operator £^ is defined as {u is set to zero) 



(10) 



7 — Tjl) 



1 



(11) 
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Starting from an initial condition {ipo{x), (po{x)) at 2; = 0, we can use Eq. (fTOl) to propagate 
the initial condition along z. Formally, the solution may be written as 



exp{C^z) 



00 



(12) 



Apparently 7 is a eigenvalue if and only if there exists an initial condition at 2; = such 
that when propagated to both end caps by Eq. (fT2l) . the boundary conditions are satisfied. 
In the present case, we may consider the linear mapping 





1 — > 


00 





^[X, z 

fx, z = 



L/2) 
-L/2) 



(13) 



To simplify the notation, hereafter we will use xo to denote the initial condition. To satisfy 
the boundary conditions, we require 

I^7,LXo| 



a^^L = mm 



IXol 



0, 



(14) 



where 1 1 . 1 1 is some suitable norm. Here we use the straightforward generalization of the usual 
norm: 



IXol 



^0 



\Si,lXo\ 



dx. 



{x, -L/2) 



dx. 



(15) 
(16) 



0(x,L/2) H 

Now if the system length L is given, one can scan a^^i over the whole space of complex 7. 
Those 7's where a^^i = are the eigenvalues. Alternatively one may fix 7 and scan through 
L to find out the system lengths for which the given 7 is an eigenvalue. 

To implement this idea numerically, first we need an accurate discrete app roximation for 
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the operator C^. Pseudospectral methods are suitable for the task. [38 
integrations for Eqs. ( |T5l) and ( |T6i) are done with Gaussian quadrature. [41 
dealing with an infinite domain as most textbooks do, we need to map (—00, 00 
to (—1,1) to apply pseudospectral methods. We define the mapping as follows. 



|The numerical 



421 1 Since we are 
smoothly 



X 



CiXi 



X, 



maps Xi G (—1, 1) to x G {—00, 00). Then a second mapping 



X2 = (1 - cs)xi + C3- tan ^ (tan (71x1/2) /C2) 



43| First 
(17) 

(18) 
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maps Xi G (—1,1) to X2 G (—1,1). This mapping can further concentrate resolution around 
X = to resolve the resistive layer. The constants ci, C2, C3 are free parameters. A set 
of Legendre-Gauss collocation points 401] are set up in X2- The differentiation matrices 
and Gaussian quadrature weights are first calculated in X2 with respect to these nodes, then 
transformed back to x. In this discrete approximation, the initial condition x is represented 
by a column vector g, which consists of 0o and ipo on the collocation points. The operators 
and Q.^^L are represented by matrices and G^^l, respectively. The norm is represented 

by 

q^Wq, (19) 



I I |2 

m 



where f denotes the conjugate transpose and is a real diagonal matrix from the Gaussian 
quadrature. We have to find 

• {<iKLWG,,,q)y' 
^-^ = T — i^ww' — (2°) 

and the corresponding q that gives the minimum. Let W = M^M and s = Mq. Then 

{s^M-^^G\ ^MWG^ lM-^sY/^ 

^7,^ = ' (,t,)l/2 ' • (21) 



Now we calculate the sing ular value decomposition (SVD) [4^ MG^^lM"^ = USV^, where 
U and V are unitary matrices of left and right singular vectors, and is a real nonnegative 
diagonal matrix of singular values. Then cx^ ^ is the smallest singular value, and q = M~^s, 
where s is the right singular vector (a column of V) corresponding to the smallest singular 
value. The eigenfunction can be reconstructed by propagating q along z by exp(£^z). 

Fig. [T] shows log(cr^^i) within a domain of complex 7 for the case 77 = 10~^, L = 1000. 
Those 7's where a^^i goes to zero are eigenvalues. In this particular case, we found six 
eigenvalues within the domain, all on the real axis. It is rather computationally expensive 
to scan through a complex domain. Fortunately, we found that the most unstable mode 
always has a real eigenvalue. Therefore we only scan through the real axis for most of our 
calculation. 

Here we briefly compare our method with the method of Evstatiev et al. in Ref. jl^ . 
These authors let ^2 — > A (they use ik instead of A) in Eq. ( fTOi) then solve it as an auxiliary 
eigenvalue problem. Let AjCj = £^ej, where Aj is the eigenvalue and Cj the eigenfunction of 
C^. The general solution is written as x = X]r=i ^i^^^^^i (truncated to a finite number of n 
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modes), which satisfies Eq. (fTOj) . To impose the boundary conditions, they consider 



where we use Xj to symbohcally denote the value of the solution (in the present case, the (j) 
part only) on a set of m points (which can be arbitrarily chosen, with m > n) on the end 
plates. Scanning through 7 for a given L (or through L for a given 7), the local minimum 
of a^^L corresponds to an eigenmode of the original eigenvalue problem. Now it is clear that 
the two approaches are rather similar, and the differences are only minor details. If the 
eigenmodes of form a complete basis, projecting onto the eigenbasis is a convenient and 
efficient way of finding exp(£^2;). Explicitly, if Xo = J2i ^i^iy then exp (C-yz) Xo = Xlj OjC'^'^ei. 
Instead of representing xo by its value on collocation points, one may use the expansion 
coefficients {oj}. The linear mapping Q-y^i, the norm ||.||, and a^^i can be similarly defined 
in terms of {aj}. Formulated in this way, our approach is nearly identical to theirs. The 
only real difference between between the two is that we use a^^L instead of ct^,l. It is worth 
pointing out, however, that even if the eigenmodes of are not complete, exp (C^z) may 
still exist. ^ Therefore our method does not rely on the existence of a complete eigenbasis. 

The advantage of a^^L is that it has a precise, analytical definition, Eq. (fT4|) : and Eq. (!20|l 
is merely a numerical approximation to it. This allows us an easier way to check convergence 
— if the numerical solution is fully converged, a^^i should be independent of the resolution 
and the mapping parameters. On the other hand, the definition of a^^i depends on the points 
chosen to enforce the boundary condition, the set of eigenmodes {ei}^^i for expansion, 
as well as the normalization of the eigenmodes. Fig. [2] shows a case of our convergence 
test with respect to the resolution, for = 50, 75, 100. The = 50 calculation gives 
7 ^ 1.301 X 10~^, while the other two calculations give 7 ^ 1.299 x 10""^. Note that the 
curves from N = 75 and = 100 calculations are almost identical, indicating convergence. 
For > 100 calculations, the results are virtually indistinguishable from the A^ = 100 
curve. Fig. [3] shows the eigenfunction ip from the same calculations. Although the relative 
error in the eigenvalue from the A^ = 50 calculation is only about 0.2%, the eigenfunction 

^ As a simple example, consider A = [3, —2; 2, —1], which has only one eigenvalue 1 as a double root and 
only one corresponding eigenvector (1,1); therefore A does not have a complete eigenbasis. However, 
exp(Az) = [1 + 2z, — 2z; 2z, 1 — 2z] is well defined. 




(Er=i h 



(22) 
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is clearly not accurate. The N = 75 result is much better, but there are still wiggles on the 
imaginary part. The = 100 calculation gives a fully converged eigenfunction. This once 
again demonstrates that an accurate eigenfunction is much more difficult to obtain than the 
eigenvalue. It took significant effort to make sure all the numerical results presented in this 
work are fully converged. 



V. RESULTS 

We use the Harris sheet By = 6tanh(x), b = 0.1, as the equilibrium. The wavenumber 
ky = 0.5 is used for most of the calculations. The dependence on ky will be addressed later. 



A. Periodic System 

The periodic system is usually solved for the case k^ = 0, but the analysis can be easily 
generalized to /c^ 7^ cases. In the so-called constant approximation (i.e. ijj is approxi 
mately constant within the resistive layer), the growthrate 7 is proportional to r]^^^. 36l. |37 
Fig. m shows the growthrate 7 divided by r/^/^ as a function of k^, for different 77. In the 
asymptotic limit 77 — 0, all curves should lie on top of each other. Some interesting features 
may be readily observed: (1) Except for the case 77 = 10~^, /c^ = is not the peak of the 
growthrate. As we will see, this has an interesting consequence in line-tied solutions. (2) 
The 7 ~ ?7^/^ scaling works better around k^ = 0. This may be due to the fact that it is 
easier to satisfy the constant ip approximation there. However, we expect that as long as rj 
is sufficiently small, eventually the constant ip approximation will work for all kz- 

At first sight, it may seem counter-intuitive that fc^ = is a local minimum of 7. In the 
tearing mode analysis 7 is determined by A' and (k ■ B)' at the resonant surface. For the 
Harris sheet profile A' may be difficult to calculate analytically when kz 7^ 0. However, for a 
similar piecewise linear profile, i.e. By = bx when < 1 and By = bx/\x\ when > 1, an 
analytic A' can be found. It may be shown that in this case fc^ = is also a local minimum. 
In fact. A' diverges as kz ±kyb. This pathological behavior may be due to the abrupt 
turning of By at x = ±1, which is unphysical. However, the result from this simple system 
suggests that it may be generic, or at least rather unsurprising, that = is the local 
minimum for the system we consider here. 
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B. Line-tied System 



The growthrate 7 of the fastest growing mode as a function of L for various 77 is sum- 
marized in Fig. O In the plot we normahzed 7 with respect to 77. As discussed in Sec. IIIH 
we expect the critical length Lc to be independent of rj, and 7 to be proportional to rj near 
marginality. Both are borne out by the numerical calculations. It is found that Lc = 115.09 
for all r], and the 7/77 versus L curves all coincide with the prediction from the force-free 
approximation. As L becomes longer, the curve starts to deviate from the force-free approx- 
imation. The deviation occurs earlier the smaller 77 is. This is also what we expect in Sec. 
mil The horizontal bars denote the range determined by the periodic 7 at /c^ = and the 
maximum periodic 7 (see Fig. Hj). We observe that as the growthrate peels away from the 
force-free approximation, it first goes beyond the periodic 7 at = 0, and eventually levels 
off at (but never exceeds) the maximum periodic 7. Therefore the force-free approximation 
and the maximum periodic growthrate may be regarded as two asymptotes of the line-tied 
growthrate. We are not able to achieve fully converged results as the t] = 10"^ and 10~^ 
cases approach their periodic limits. However, based on the observation for larger 7/ cases, 
it is fairly safe to say that the two-asymptote scheme may be generally applicable. 



and 



for various L, with 7/ = 10 



Fig. [6] and Fig. [7] show the eigenfunctions 
The eigenfunction is normalized such that ^/;(0, 0) = 1. The eigenfunction is very smooth at 
the critical length L^. = 115.09. As L becomes longer, the internal layer at x = become 
steeper and steeper. Up to L = 1000, the eigenfunction still largely covers the whole length 
along z. At L = 1500 and 2000, the eigenfunction becomes localized to the midplane. This 
localization is more pronounced in ip. At L = 3500, the growthrate is greater than the 
periodic one at kz = 0. The internal layer "splits" into two, and there are wiggles along the 
z direction. This is due to the peculiar two-bump structure in the dispersion relation of the 
periodic system (Fig. Hj). In the periodic system, 7 > 'j\k^=o corresponds to four nonzero 
fc^'s, therefore we expect the line-tied mode to be dominated by ^2 7^ behavior. The four 
periodic modes have resonant surfaces at four different locations. The four resonant surfaces 
form two pairs, one in x > and the other in x < 0. These two pairs of resonant surfaces 
roughly correspond to the two parts of the split internal layer. From the discussion in Ref. 
jl^. I15I. [igI , we expect an expansion from four eigenmodes of to be a good approximation 
when the two resonant surfaces within each pair become sufficiently close to each other. 
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However, it is found that for this to occur one would have to go to an enormously large L, 
and we have trouble achieving full convergence in this regime. 

As a comparison, Fig. [S] shows the eigenfunction with the same parameters as Fig. 
[HI but in the force-free approximation. The force-free approximation for L = 400 is almost 
identical with the real solution. Even at L = 1000, the force-free solution is still not too far 
from the real one. At L = 1500 and beyond, the difference between the two is significant. 
The same conclusion can be drawn from Fig. [5l where the force-free dispersion relation is 
good for up to L = 1000, then it starts to depart from the real one. We also observe that the 
internal layer of the force-free eigenfunction keep getting steeper as L increases. Fig. [9] shows 
the midplane cut of lm(0) for various L, and Fig. [TU] shows its force-free counterpart. Also 
shown in both figures is the kz = periodic eigenmode with the inertia term included, for 
reference. Notice that in periodic systems, the inertia term has to be included to resolve the 
inner layer, 36|, |37(| while in line-tied configurations the inner layer can be resolved with the 
help of the boundary condition (no k ■ B = singularity) alone. Following the convention of 
Ref . [l^ , the internal layer width of the force- free approximation is called the "geometric" 
width, and that of the periodic eigenmode is called the "tearing" width. At L = 1500, the 
geometric width is approximately the same as the tearing width, and this is exactly when 
the force-free approximation is about to fail. 

Therefore, we come to the conclusion that the transition from the resistive scaling 7 ~ 77 
to the tearing scaling 7 ~ r]^^^ occurs when the inertia term is no longer negligible, or 
equivalently, when the tearing width is comparable to the geometric width. This prompts 



us to do a scaling analysis similar to the one of Ref. 16|]. Assuming the constant ijj 
approximation, and letting S be the thickness of the internal layer, we have V(f) ~ 0/^^, 
Vtp ~ ^AV(5, and By ~ B' 6 within the layer, where A' is the jump of ip' across the layer 



divided by '0.[36|,l37| For the periodic case, we have ■yVcp ~ kyByVip and '-yip ~ kyByCp ~ rjVijj 



within the layer. That is 

70/5' ~ kyB'^A'iP, (23) 
77/i ~ kyB'y6^ ~ 7]tljA'/S. (24) 
From that, we find the tearing growthrate 

7tear^ns ~ {kyB'^/^A'^^'fj'/', (25) 
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and the tearing width 

^tearing ~ 7^"" {kyB'^r^'K (26) 

For the hne-tied case with the force-free approximation, we may assume dz ^ ^/ L since the 
eigenfunction spreads over the entire length (Fig. [H]). Within the internal layer, we have 
dz ~ kyBy] i.e., the geometric width 

^geometric ^ 

1/kyB'^L. (27) 

From the induction equation, we have 

Iforce-free ~ V^'/^ ~ r]kyB'yLA' . (28) 

Also, we may estimate from ■yip ~ kyB'^Scp as 

~ 7]kyB'yL^A'ip. (29) 

The force-free approximation is self-consistent if yVcj) <^ dz'Dip; i.e., ycp/S'^ ^ tpA'/SL. 
Using Eqs. fl27j) . fl28l) . and fl29l) . the self-consistency condition can be written as 

L « Li ~ r/"2/5A'-i/5(A;,5;)-3/5, (30) 

where is the transition length from resistive scaling to tearing scaling. The same condition 
can be obtained if we require Stearing < ^geometric- 

The scaling analysis predicts that the transition length Lt ~ r/"^/^. Also we know that 
as L becomes large, 7 ~ 77^/^. This motivates us to rescale Fig. [5] in the following way. We 
plot 7/77^/^ versus (L — Lc)'rf^^. Normalizing 7 by r]^^^ brings the asymptotic regime to the 
same level, while multiplying (L — Lc) by 77^/^ brings the transition length to the same place. 
The result is shown in Fig. [TTl Clearly those curves for different rj become close to each 
other, especially those with smaller 77 are almost identical. Therefore the numerical results 
agree with the Lt ~ 77-2/5 prediction. The analysis also predicts y force- free ~ rjkyB'yLA' . 
We can see from Fig [5] that the force-free dispersion relation is very close to a straight 
line. That suggests that we may associate the slope with y force- free/ ~ ^yB'yA' . For the 
Harris sheet profile. A' = 2(1/A;,y — ky). To test the conjecture we have to find the slope for 
different ky. The slope and the critical length for different ky are summarized in Table [Til 
We found that Lc only weakly depends on ky. Fitting the slope with const x kyA'B'y gives 
slope ~ 0.0159kyA'B'y. This turns out to be an excellent fit, as shown in Fig. [T2j 
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ky Lc Slope kyA'B'y 
0.3 112.85 0.0288 0.182 
0.4 111.50 0.0266 0.168 
0.5 115.09 0.0239 0.150 
0.6 123.86 0.0204 0.128 
0.7 141.09 0.0162 0.102 

Table II: The critical length L^. and the slope of the force- free dispersion relation 7/77 versus L 
for different ky. Also shown is the corresponding kyA'B'y. The scaling analysis suggests that slope 
~ kyA'B'y. 

VI. SUMMARY AND DISCUSSION 

In this work, we study the effects of line-tying on the resistive tearing instabihty in great 
detail. We found that hne-tying has a stabilizing effect, and instability occurs only when 
the system length is greater than a critical length which is independent of the resistivity 
77. When the system length is not too long, the plasma inertia is negligible. In other words, 
the plasma is approximately force-free throughout the whole domain. In this regime, the 
growthrate 7 is proportional to rj and the internal layer of the unstable mode is much wider 
compared to the periodic mode. The tearing mode scaling 7 ~ 77^/^ is recovered when the 
system becomes sufficiently long. The transition from 7 ~ 77 to 7 ~ 77^/^ occurs when 
the plasma inertia is no longer negligible within the resistive layer; this is also when the 
"geometric" width becomes comparable to the "tearing" width. The transition length Lt 
scales like rj'"^^^. 

Physically it is not hard to understand why the growthrate in the force-free approximation 
is proportional to 7/. In the original system, there are two time scales, the Alfvenic time scale 
and the resistive time scale; and the tearing mode grows at a hybrid time scale of the two. If 
the inertia is neglected completely, the resistive time scale is the only one left; therefore the 
mode can only grow on the resistive time scale. However, neglecting the inertia in the entire 
domain is possible only in a line-tied system, as the line-tied boundary conditions serve as a 
means of resolving the singularity at the rational surface. In a periodic system, the inertia 
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term always comes into play in the vicinity of the rational surface. 

The periodic growthrate and the growthrate from the force-free approximation may be 
regarded as two asymptotes of the line-tied growthrate. The force-free growthrate 'j force-free 
is approximately linear with respect to L and the slope is found to be related to A'. To 
a very good approximation force- free — -0159 ky A' Byr]{L — L^), where B'^ is evaluated at 
X = 0. Unfortunately, we do not have a good way to estimate the critical length Lc- 

To a large extent, our findings are in agreement with the earlier study by Delzanno and 
Finn in cylindrical geometry. 15|] One major discrepancy is that we found Lt ~ 7^~2/5 ^j^jjg 
they reported Lt ~ r^^^/'^-^ ^-0.24 (xhey actually consider the crossover resistivity rjcross 
for a given L, and find rjcross ~ L"^'^; but this is just a matter of which variable is held 
fixed). A possible explanation for the discrepancy is that they use a fairly large viscosity, 
which remains unchanged when other parameters are varied. We have not studied the effect 
of viscosity in great detail. Some preliminary results indicate that viscosity does have a 
significant effect, and Lt ~ 7^"^/^ scaling is no longer valid when a large viscosity is added. 

We conclude with a brief discussion of our findings in the context of solar corona. The 
solar corona is very highly conducting, and the aspect ratios of observed coronal loops are 10 
-100, far smaller than necessary for growthrates in the tearing regime. Therefore it is likely 
the resistive instability will operate in the 7 ~ 77 regime. An immediate question is that, 
since the resistive time scale is so long in the solar corona, is the resistive instability relevant 
at all? We do not have a clear answer as yet. It is quite clear that our model is highly 
idealized and has some obvious shortcomings. An immedaite concern is that footpoint 
drivers are not included in our analysis; therefore the equilibrium profile will decay on a 
resistive time scale r^. Our analysis implicitly assumes either that 7 ^ 1/Tr?, such that the 
resistive decay may be neglected, or that the equilibrium is maintained by some external 
current drive, as in laboratory devices. In the solar corona the electric current is driven 
by footpoint motions, and the assumption is clearly inadequate when Lc < L Lt, as the 
growthrate 7 is either comparable to or even smaller than l/r^. A remedy to this problem 
would be adding a footpoint driver to the system to maintain the equilibrium. This will 
add an overall background flow to the equilibrium and potentially have some effect on the 
stability, especially when the system is close to marginal stability. Detailed analysis is left to 
a future study. Furthermore, the rigid wall boundary condition, which we take for simplicity, 
may be too restrictive as argued by some workers. This is especially true when considering 
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modes with 7 ~ r], as the photosphere may be able to respond during the slow time scale of 
growth. A proper treatment requires some detailed knowledge of the photosphere interior. 
7, 4^, 45, 46] In light of these shortcomings, our findings certainly need furthur scrutiny 



when applied to the solar corona. 

As a final remark, the nonlinear evolution of the mode has to be studied. Indeed, even 
in the classical tearing mode theory without line-tying, the unstable mode evolves into the 



nonlinear regime very quickly. 



471 ] Our preliminary study indicates that thin current filaments 



can be created within the original current layer as the mode evolves nonlinearly. Therefore 
the role resistive instabilities play in current sheet thinning cannot be overlooked. Another 
interesting point is that, since the stability criterion does not depend on the resistivity but 
the growthrate does, the solar corona may appear to be quiescent but already be in an 
unstable configuration. The slow growthrate of the resistive instability allows the corona to 
be driven well beyond the critical point, therefore more free energy could be stored. Yet the 
resistive instability, however slow it may be, also provides a means for the corona to tap into 
the free energy in an otherwise ideally stable configuration. The multiple roles the resistive 
instability may play here leave much to ponder about. We hope to address some of these 
issues in the future. 
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Figure 1: (Color online) log{a^^L) in a domain of complex 7. = 10"^ and L = 1000. Those 7's 
where a^^i goes to zero are eigenvalues. Here we found six eigenvalues within this domain. All of 
them are on the real axis. 
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N=100 



1 .295 1 .3 1 .305 

T X10-' 

Figure 2: (Color online) Convergence test for the case 77 = 3 x 10~^, L = 2000. The mapping 
parameters are ci = 2, C2 = 0.1, C3 = 0.9. The = 50 calculation gives 7 ~ 1.301 X 10""^, while 
the other two calculations give 7 ~ 1.299 X 10"^. Note that the curves from N = 75 and = 100 
calculations are almost identical, that indicates convergence. For > 100 calculations, the results 
are virtually indistinguishable from the N = 100 curve. 
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Figure 3: (Color online) The real and imaginary parts of eigenfunctions from the convergence 
test shown in Fig. [2j (a) N=50, (b) N=75, (c) N=100. Note that although the relative error on 
eigenvalue from = 50 calculation is only about 0.2%, the eigenfunction is clearly not accurate. 
The N = 75 result is much better, but there are still wiggles on the imaginary part. The iV = 100 
calculation gives a fully converged eigenfunction. 
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Figure 4: (Color online) as a function of k^, for periodic systems. 
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Figure 5: (Color online) 7/77 as a function of L, for various r]. The horizontal bars denote the range 
determined by the periodic 7 at /c^ = and the the maximum periodic 7 (see Fig. [3]). The dashed 
line is the prediction from the force-free approximation. Inset: the expanded view of the lower left 
corner. 
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Figure 7: (Color online) The eigenfunction for various L. t] = 10 
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Figure 9: (Color online) The midplane cut of lm(0). Also shown is the periodic solution with 
kz = for reference. 
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Figure 10: (Color online) The midplane cut of lm(0) for the force-free approximation. Also shown 
is the kz = periodic solution (with inertia) for reference. 
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Figure 11: (Color online) Rescaled line-tied growthrate as a function of L. 
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Figure 12: Slopes of the force-free 7/77 versus L curve, which is approximately a straight line, for 
different ky. The curve 0.0159/cyA'i?y is an excellent fit with the numerical results. 
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